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Abstract 

The field theory of a short range spin glass with Gaussian random interactions, is 
considered near the upper critical dimension six. In the glassy phase, replica sym- 
metry breaking is accompanied with massless Goldstone modes, generated by the 
breaking of reparametrization invariance of a Parisi type solution. Twisted boundary 
conditions are thus imposed at two opposite ends of the system in order to study the 
size dependence of the twist free energy. A loop-expansion is performed to first order 
around a twisted background. It is found, as expected but it is non trivial, that the 
theory does renormalize around such backgrounds, as well as for the bulk. However 
two main differences appear, in comparison with simple ferromagnetic transitions : 
(i) the loop expansion yields a (negative) anomaly in the size dependence of the free 
energy, thereby lifting the lower critical dimension to a value greater than two given 
hj dc = 2 — ri{dc) (ii) the free energy is lowered by twisting the boundary conditions. 
This sign may refiect a spontaneous spatial non-uniformity of the order parameter. 
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1 Introduction 



Spontaneously broken symmetries are characterized by the existence of several pos- 
sible pure states. If one imposes "twisted" boundary conditions , i.e. different pure 
states at two ends of the system, the free energy per unit volume will be slightly 
greater than the free energy corresponding to one single pure state over the whole 
system. 

For a simple discrete symmetry, such as the Z2-symmetry of Ising-like systems, one 
may consider an (hyper)-cubic system with up spins in the z = plane, down spins 
in the z = L plane and for instance periodic boundary conditions in the transverse 
directions Xi,X2,- ■ - xa-i- This will generate an interface in the system centered 
around some plane z = Zq and a cost in free energy 

AF = - = aL"-' (1) 

in which cr(T) is the interfacial tension. As is well-known the power {d — 1) of L in (|ip 
implies that the lower critical dimension of systems with a discrete symmetry is equal 
to one, i.e. there is no ordered phase unless d is greater than one. At leading order 
the classical (mean field) configuration for the order parameter, given the boundary 
conditions, is a kink of hyperbolic tangent shape, interpolating between up and down 
spins. The fluctuations are given at one-loop order by the Fredholm determinant 
of a one-dimensional Schrodinger operator in a l/cosh^(2; — Zq) potential which, 
as is well-known, is solvable analytically. Every term in the loop expansion for the 
free energy about mean field theory, is then proportional to L'^^^, and the succesive 
contributions build up the correct exponent and amplitude for the interfacial tension 
cr. 

For continuum spontaneously broken symmetries, the situation about the upper 
critical dimension is technically different. For an N- vector model one considers for 
definiteness an order parameter, which is uniform along the vector (1, 0, ■ ■ ■ , 0) in the 
z = plane, and uniform but rotated by an angle 9o in the plane z = L, i.e. lying 
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along the vector (cos 6'o, sin 6'o, 0, ■ ■ ■ , 0). There again one expects a cost in free energy 

AF = (T{T,eo)L'^-^ (2) 

in agreement with a lower critical dimension equal to two, and with a "twist" energy 
(t(T, ^o) (or spin stiffness constant) vanishing as 6'q for small ^o, ( the ratio ct/Oq is 
the helicity modulus |^). If it is quite elementary to verify these statements within 
mean field theory, not difficult also to check them in the vicinity of the lower critical 
dimension di = 2 through the non-linear sigma model ^. Near the upper critical 
dimension du = 4, things are not as simple. The mean field solution is not elementary 
and one may fear that the loop expansion might be difficult to handle. However it 
turns out that for L large, the analysis of fluctuations is simply perturbative and 
finally explicit. It follows from this analysis that the massless Goldstone modes give, 
as expected, an L''"^ behaviour in the twist free energy to all orders in the loop 
expansion. 

For a spin glass the nature of the broken symmetry in the low temperature phase 
is more difficult to visualize. However within the replica approach, and Parisi's 
ansatz for the mean field solution there are indeed "replicon" massless Goldstone 
modes ^ (plus "anomalous" massless modes). The broken symmetry at the origin of 
those modes may be related to a reparametrization invariance of the action. More 
specifically the mean field solution depends, in the continuum limit of Parisi's scheme 
of replica symmetry breaking for the Edwards-Anderson model, of two functions p{t) 
and Q{t) in which t is the continuum labelling of the steps of breaking, Q{t) the Parisi 
order parameter and p{t), the continuum limit of the size of the successive boxes in 
which the n replicas are divided. The free energy is not a separate function of Q{t) 
and p{t) but depends only of Q as a function of p, leading for instance to the simple 
"gauge choice" p{t) = t of Parisi. The existence of massless modes may be related to 
this arbitrariness ||l2|. However one does not see any physical "external field" which 
could be used to tune a given specific gauge choice. The situation is thus reminiscent 
of cases such as superfluid Helium, in which there is no physical conjugate variable 
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to the order parameter which one could use to fix its phase. However if one takes 
two samples, they have no reason to carry the same phase, and this phase difference 
manifest itself in Josephson's junctions for instance. 

In this note we report the result of an analysis, in which one imposes again two 
different schemes at two ends of the system. In the z = plane we have chosen the 
simple Parisi gauge 

p{t,z = 0)=t (3) 
whereas in the z = L plane we have imposed 

p{t, z = L)=t + hit) (4) 

in which we assume that h{t) is some given infinitesimal function, vanishing with t, 
with support < t < x. All calculations have been performed to lowest order in 
h{t). The mean field solution, to lowest order in h{t), provides a linear interpolation 
beteween the two end planes, and a free energy which is proportional to L"'"^ as for 
the N-vector model. At one-loop order, in dimension d = 6 — e one finds after a 
long calculation, whose details will be reported elsewhere, a free energy for the twist 
which is proportional to L'^"^ logL. Those logarithms, which are caused here by 
the absence of a mass gap to the Goldstone modes, change drastically the situation 
compared with ordered states. They may be exponentiated in the standard way and 
yield, to first order in e a twist free energy which is proportional to 

Jo 

in which r measures the temperature below the glassy transition. This is, up to 
one-loop, the approximation to 

Since rj is negative, this shows that the lower critical dimension dc is larger than two, 
given by 

4 = 2-r/(4) (7) 
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The sign of the result is a puzzle on which we have a few comments at the end. 



2 Mean field theory 

The action for the Edwards- Anderson spin glass is written in terms of an n x n matrix 
Qab, in which a and b are replica indices and 

In a Parisi replica infinite symmetry breaking scheme, one divides the n — po 
replicas into po/pi boxes of size pi, each box of size pi is divided into pi/p2 boxes of 
size p2 , and so on, ad infinitum. The matrix elements Qab follow those steps and are 
characterized by a correlative infinite sequence Qo,Qi, - ■ ■ ■ 

In the continuum limit, wc arc thus led to an action which depends on two spatially 
varying functions p{t,z) and Q{t,z), in which < z < L and t refers to the steps 
in the symmetry breaking scheme. ( In the {d — 1) transverse directions, periodic 
boundary conditions have been imposed, and the mean field solution is independent 
of those tranverse space variables) . In terms of those functions the action reads 

- ^ dtp{t, z){p(t, z)Q\t, z) + 3g2(t, dsp{s, z)Q{s, z)'^] (9) 

( p, Q denote derivatives with respect to t) . 

In the bulk, with non twisted boundary conditions, this action is manifestly a 
function of Q{p) alone, and not separately of Q{t) and p{t). Indeed the extrema of 
this free energy are given as solutions of 

A{t) =tQ + ^Q'-^ (pQ' + f ds^Q'is) + 2Q [' ds^Q{s)] = 0. (10) 
6 2 \ Jo ds Jt ds J 

and Parisi 's solution is 

w 

Q{t) = ^^'(^) 0<t<xi 

Q{t) ^ Qi for xi<t<l (11) 
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with the Edwards-Anderson order parameter Qi defined by 

T + uQl-wQi=0. (12) 

With the twisted boundary conditions we obtain two variational equations for 
Q{t, z) and p(t, z), which read 

1 d d'^O 

M+4^(pQ) + ^])^ = (13) 

in which A is defined in (p!0|). 

To lowest order in the imposed twist h{t) one checks easily that the solution is 

2u z 
p{t,z) = —Q{t,z)=t+-h{t) (14) 
w L 

for < t < x (x is the end of the support of h{t) and then the solution is the bulk 
one for t > x. To second order in h{t) though, the bulk proportionnality of p and Q 
is lost). 

The incremental free energy, which follows from this twisted solution (compared to 
the bulk one), comes purely from the kinetic energy (since the bulk relation p(t, z) = 
^Q(t,z) still holds). The final mean field result is, at lowest order in h{t), 

1 /in rti 

AF,_, = -lim-(Z,';,3,,,-ZbV) = -i-fL"-' dth{t){h{t)+th{t)) 
n^o n 2u Jo 

This result holds above the upper critical dimension du = Q but, contrary to ordered 
states, we shall see that it is modified by fluctuations. We now proceed to the one-loop 
computation. 

3 One-loop fluctuations around mean fleld 

The theory is now extended in dimension d = 6 — e, dimensionally regularized , and 
later renormalized. 
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1. Consider first the Replicon sector where the fiuctuation matrix may be 
fully diagonalized ( for a review of fiuctuations beyond mean field see 
[§]). Its continuation to the bulk would write, in the continuum limit, 



p(R) 

loop 



1/ 

2 



1 dt dk d r^dl d 
'2 Jt TdkJt Idl 



d'^p 



where we have used the notation 
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w 
2u' 



(16) 



(17) 



In (|T6|) the first bracket comes from the multiplicity of the replicon 
modes with associated eigenvalues p"^ + Ao{k, I; t) as in the argument 
of the logarithm 

Ao(fc,/;t) = |(fc^ + /2-2t^). (18) 

Under the twist (p!^, the above bulk result (|16|) is changed in two ways. 
First the argument of the logarithm is to be replaced by 



(19) 



with 



Ai = g[kh{k) + lh{l) - 2th{t)] (20) 

In fact there is also a quadratic term {z/L)^ A2, that has been omitted 
here since at one-loop it cancels through the bulk subtraction. 

Expanding now the logarithm to second order in Ai (in order to collect 
the quadratic terms in h{t)), we obtain as twist contribution the term 
L'^-^ dt /•! dk d r^dl d r d'^'^qr 



2 Jt k dkJt I dlJ {271^-^ 



d ^ f^d'^' 1+"^ ^^smKzsmKz' 



Al (21) 



In this expression boundary conditions at the two end planes z = and 
z = L have been taken into account ; indeed the fluctuating part of 
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the field vanishes at those boundaries, leading to the appropriate basis 
sin with m = 1, 2, ■ ■ -. 

Secondly, the twist ( p!^ changes also the multiplicity itself. This is 
taken into account via an identity expressing the reparametrization 
invariance under a z-independent shift Aq — Aq + Ai. As a result the 

total twist contribution is then obtained by replacing in ( pT]) ^ by 

zz' ~l/2{z^ + z'^) 

Z2 ■ 



Performing the K- integration on the modified (21) one obtains 



in which 



.1 .1 Q dl d 
~8~Vo ~2 Jt ~kdkJt Tm 

L rL p-M\z-z'\ _ -M(z+z') 
dz / dz'[- (22) 

io ^ 2M J 1 V ; 



M = ql + Ao. (23) 



Performing the z, z' and finally qt integrations, to gather poles in e and 
logarithms, one obtains : 

Ai^tSSL = -To^'-'Sa g' r dt h\t) [- + log L + ■ ■ ■], (24) 
12 JO e 

with 

2 

^ (47r)'^/2r(d/2) ^^^^ 
It is a strong argument in favour of the consistency of the calculation 
to see that the one-loop contribution is, like mean-field, proportional 
to the integral / dt h'^{t). Indeed, otherwise the fiuctuations would not 
be renormalized by a simple change in the coupling constant. In the 
intermediate steps this final form is far form obvious. In particular it 
involves the unexpected identity 

f'^ dt dk d dl d ^ n ^ . 9 , , 9. x , . 



2. In the Longitudinal- Anomalous (L-A) sector, the fluctuation matrix can 
only be diagonalized by blocks (with blocks of size (i? + 1) x (i?+ 1), 
R being equal to the number of steps of replica symmetry breaking). 
In the Parisi limit, in which R goes to infinity, the L-A contribution to 
the bulk free energy writes 

where, again, the first bracket comes from the new multiplicity of the 
L-A modes . Besides we have the (t, s) matrix 



Bk{t, s) = g(lnf{t, s))[Q{k - s) + k6{k - s) + 2e(s - k)]ds (28) 



and, as in (118 



Ao = g{k'^-f)/2 t<k 

= t>k (29) 

The calculation is more involved here, but we proceed as in the repli- 
con sector, collecting quadratic terms in h when performing the twist 
transform as in (11^) or in Aq ^ Aq + {z/L)Ai. As above we keep 
only terms quadratic in the propagator {q^ — d'^/dz'^ + Aq)"^, higher 
order terms being, at one loop, ultra-violet convergent. In contrast to 
the replicon sector (where only the twist of Aq into Aq + (z/L) Ai 
contributed) we need here to take care of the twists over the matrix 
elements Bk{t,s) and over the multiplicity. Altogether, complicated 
expressions rearrange themselves to give 

A^ii^ = \l'^-'S, / r dt h\t) [- + logL + ■ ■ ■]. (30) 
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Notice that, at one-loop, there is no contribution in log(r) i.e. involving 
the mass wQi ~ r (as in (|l^)). The reason is that neither Qi (i.e. xi), 
nor po = n, fluctuate under reparametrization. Even as a correction 
term to the contribution in (zh/L), (z'h/L), it would take to expand 
the logarithm to third order before Qi showing up. 

4 Renormalization and scaling 

Before proceeding, one has to take into account the fact that the 
^J^ab'Ptb coupling is irrelevant, and "dangerous" since the fluctua- 
tions make it singular below dimension eight 0, ||. Indeed in a pure 
^tr(j)^ theory, the one-loop contribution has the effect of replacing u by 

u^u + 12w' I , =u + 6W(2r)-(i+^/2). (31) 



We thus have for the twist free energy 
Sd 



twist 



^d-2^2+e r dth\t) \l--w\-+\ogL) + ---] (32) 

Jo 2) e 



288^6" 

( a factor Sd has been included in w^). It is now crucial to verify that the 
replacement of the coupling constant and temperature by their renor- 
malized counterpart wr,tr rids us of the 1/e poles. The computation 
of those renormalizations is easily done in the paramagnetic phase. It 
gives |TT 

TR = r[l-^ + 0(w^)]Z 
Z = [l + ^ + 0{w')] (33) 



^R 



from which follows 
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The 1/e pole is thus exactly cancelled and we end up with 

^Ftw.st = -7^L'-'r'^^^ fdt h\t) [1 - \wl\ogL + ■■■]. (35) 

AooWj^ Jo O 

If we substitute to wr the fixed point w* , zero of the /3- function 

P(wr) = -'-wr + wI + 0{w%) (36) 
one obtains that, to this order, the result exponentiates to 

288u;| ^ Jo ^ ^ ^ ^ 

Introducing the usual critical exponents rj and z/, whose e-expansions 
are known to be 



V = -le + 0{e') ^ = ^(1 + ^^ + 0(5^)), (38) 

one may write to this order 

This is reasonable since it gives the final twist free energy as a function 
of L/C, to a power, which is the expected scaling form for the ordinary 
order-disorder transitions : 

AF,^^st ~ -ijY-'-''^. (40) 
There are two differences though with ordinary transitions. 

• First the power of L/C, is non-canonical . We thus verify, to the 
order of one-loop, an extended form of scaling, appropriate when 
the soft transverse modes are not isolated, but being at the bot- 
tom of a gapless band, they are no longer infra-red free : their 
propagator develops a (negative) anomaly rj. At the lower criti- 
cal dimension dc, the twist free energy should vanish and thus dc 
should be the solution of the equation 

4 = 2-7^(4). (41) 
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This same answer had been anticipated earher on the basis of 
scahng arguments apphed to the null overlap replicon sector |[T0|| . 
Using numerical estimates for rj |T^,one gets from ( ^T]) a value 
of dc close to 2.5, whereas a self consistent mean-field approach 
for the twisted free-energy of two copies, surprisingly yields [0 
exactly 5/2. 

The sign of this twisted free energy is negative. This will be 
discussed below. 



5 Discussion 



The calculation of fluctuations around the background of a twisted 
mean field solution is renormalizable, as expected, in spite of its non 
spatial non-uniformity, as already checked in the simple 0(N)-model 
[^]. However contrary to the simple ferromagnetic transitions, the in- 
fluence of this twist on the size dependence manifests itself by a loga- 
rithmic dependence in L at one-loop which exponentiates to a negative 
anomaly. This generates an increase of the lower critical dimension. 

However another difference with ferromagnetic transitions is the sign 
dependence of the twist on the free energy : the twisting leads here 
to a decrease in the free energy. This is not in contradiction with the 
principles of thermodynamics. A similar situation could occur with an 
antiferromagnet, since there as well, the free energy may be lowered by 
imposing external fields of opposite signs at two ends of the sample. 
In such a circumstance, one expects that the system would sponta- 
neously breaks spatial uniformity and develop a space dependence in 
the "gauge" choice of the replica symmetry breaking. Unfortunately 
we are not aware of any conjugate field to the gauge choice p{t,z) which 
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could be imposed to improve the mean field starting point. Transpos- 
ing to dynamics the instability to twisting found here, might be the 
sign of a space dependence of local time scales, as discussed in recent 



simulations of finite dimensional spin glasses . 
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